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Summary. — The leading hadronic effects in electroweak theory derive from vac- 
uum polarization which are non-perturbative hadronic contributions to the running 
of the gauge couplings, the electromagnetic aem(s) and the SU{2)l coupling 02(5). 
I will report on my recent package alphaQED [1], which besides the effective fine 
structure constant acm(s) also allows for a fairly precise calculation of the SU{2)l 
gauge coupling 02(3). I will briefly review the role, future requirements and possi- 
bilities. Applied together with the Rhad package by Harlander and Steinhauser [2], 
the package allows to calculate all SM running couplings as well as running sin^ S 
versions with state-of-the-art accuracy. 

PACS 12. 15. -y, 13. 40. -f, 12.15. Lk, 13.40. Ks - . 



1. Introduction 

Precise Standard Model (SM) predictions require to determine the U {1)y ^S) SU (2) l^ 
SU{3)c SM gauge couplings aem, ci2 and a^ = as (QCD) as accurately as possible. 
Obviously, the predictability of theory is limited by the precision of its input parameters. 
This in particular requires to fight precision limitations due to non-perturbative hadronic 
contributions. Precise predictions confronting precise measurements are the basis for all 
SM precision tests, which allow us to unravel new physics from discrepancies between 
theory and experiment. An important test case, which requires as precise as possible 
running couplings, is the quest of gauge coupling unification in grand unified extensions 
of the SM. 

Key input parameter for ILC physics currently are known to precision: 
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We observe that the accuracy of a(Af2) is roughly one order of magnitude worse than that 
of the next best Mzl The loss in precision caused by non-perturbative strong interaction 
effects is 10^ between the classical low energy a and a{Mz). The requirement for ILC 
precision physics is 

(2) ^^-5x10-. 

A prominent example where theory may be obscured by lack of precision in the effective 
a is the indirect Higgs mass bound obtained from the precise measurement of sin 0^^. 
The required improvement could be achieved by dedicated efforts in cross-section mea- 
surements in the energy range from 1.2 to 3.2 GeV, and be adopting the Adler function 
controlled split in parts evaluated from data (from experiments or from lattice QCD 
simulations) and parts which can be calculated reliably in perturbativc QCD (pQCD): 



, (5) /,,f2\ _ A^(5) / Nidata 



^<Zm) - Aa\::,i-s„r^^ + Aa- (-A/^) - Aa- (-.o) 



,,(5) I- nr2\ A„,(5) 



pQCD 



(3) Aa^^^,{Ml)-Aal^^,i-Ml) 



pQCD 



where sq can be optimized by adopting the Adler function as a monitor for the range of 
validity of pQCD [17, 18]. In the following we will present a description of the package 
alphaQED which allows state-of-the-art calculations of the SM running couplings, option- 
ally, with their imaginary parts. Some emphasis is put on the not so straight forward 
determination of the running SU{2)l coupling a2(s), which is important for the calcula- 
tion of variants of the weak mixing parameter sin Qwis), an interesting quasi-observable 
and monitor of new physics particularly at ILC energy scales. 

2. Effective running coupling oqed 

The effective fine-structure "constant" a{E) depends on the energy scale because of 
charge screening by vacuum polarization: 

(4) Aa{s) = -e^ [Re n"^''(.s) - n''''^(0)] 

which exhibit the leading hadronic non-perturbative part Aj^^^a. n(s) = 11(0) + sn'(s) 
denotes the transversal current correlator, for the electromagnetic current 11(0) ~ 0. 
While electroweak effects (leptons etc.) are calculable in perturbation theory, the calcu- 
lation of the strong interaction effects (hadrons/quarks etc.) by perturbativc QCD fails. 
Fortunately, dispersion relations and the optical theorem allow us to perform rather 
accurate evaluations in terms of experimental e"*"e~-data encoded in 

jy , ■._ (J{e+e- ^ 7* ^ hadrons) 

(5) K^{s) = — -—^ ' — . 

^' ^^ ' a{e+e- ^j* ^ fi+fi-) 
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For the electromagnetic running coupling the dispersion integral reads 

TT cut 

The high energy tail is neatly calculable perturbatively by the virtue of asymptotic 
freedom of QCD. Errors of data imply theoretical uncertainties. Some of the data sets 
are old and of rather limited precision, especially in the range above 1.4 GeV to about 
2.2 GeV, a range which is subject to new measurement at the VEPP 2000 facility at 
Novosibirsk. Data from different experiments are combined by standard methods as 
recommended by the Particle Data Group (see e.g. [3]). In recent years progress has 
been due to much better a{e~^e~ — > hadrons) determinations at Novosibirsk (CMD2, 
SND) [4, 5] and more recently by the novel radiative return high accuracy measurements 
by KLOE [6, 7] and BABAR [8] (see also [9, 10]). Typically, vacuum polarization leads to 
large corrections and in fact a{E) is steeply increasing at low E already. So the deviation 
of a(77ip) at the muon mass scale jti^i from a gives the big leading hadronic correction to 
the muon g ^ 2 [11]. That is why we need to know the running of aqED very precisely 
at all scales (see Fig. 1). Non-perturbative hadronic effects in electroweak precision 
observables affect most SM predictions via non-perturbative effects in parameter shifts, 
typically: 



(7) sin" 8i cos" 8 



V2 G^ A/| 1 - Ari 
where 

(8) Ari = Ari(a, G^, Mz^mn, ruf^t, rrit) 

represent the quantum corrections from gauge boson self-energies, vertex- and box- 
corrections. Uncertainties obscure in particular the indirect bounds on the Higgs mass ob- 
tained from electroweak precision measurements. Basic observables like Mw [sin Qw = 
1 - M^/M|], g2 [sin^Qg = e'^/g^ = {Tra)/{V2G^ M^)] or the vector coupling Vf 
[sin^ Qf = (4|Q/|)^^ (1 — Vf/ttf) , j ^ v\ are related to versions of sin^ 8vk obtained 
form (7) and the general form of Ar^ reads 



(9) Ar, = Aa - /.(sin^ 6,) Ap + Ar, , 



cinaindcr 



with a universal term Aa which affects the predictions for Mw ^ ^Lfi, ^fb^ ^ i^ ^**^- Only 
the p parameter in the axial coupling a/, which is renormalized by pf ~ 1/(1 — Ap), is 
independent from leading non-perturbative hadronic effects. 

One issue concerning running couplings concerns the question complex vs. real a{s). 
For s 7^ (4) provides the definition of a complex coupling if we relax from taking the real 
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Fig. 1. - Aacm(-E) and Aa2{E) as functions of energy E in the time-like and space-like domain. 
The smooth space-like correction (dashed line) agrees rather well with the non-resonant "back- 
ground" above the (ji-resonance (kind of duality). In resonance regions as expected "agreement" 
is observed in the mean, with huge local deviations. 



part only. A typical example where this matters is the vacuum polarization correction to 
be performed on R{s) before it can be used in (6): i?physicai -^ R = (c«/'^(s))^ ^physical- 
Usually, a{s) is take to be real, i.e., (q;/q(s))^ = |1 — Ren'(s)p (n'(0) subtracted). 
More precisely, one should subtract |1 — n'(s)p ~ a/\ac{s)\)'^ where a^s) denotes the 
complex version of running a. Typically, corrections from imaginary parts given by 
1 — |1 — n'(s)p/(Q;/a(s))^, are small ^ 0.1 % in non-resonance regions. However, at 
resonances corrections are of order ^ l/Tpt and thus are large for narrow resonances. 



3. The coupling a2, Mw and sin^ 8/ 

Unlike for the electromagnetic coupling, for the SU{2)l coupling the hadronic shift 
cannot be directly obtained by integration of measured data. There is however a pretty 
clean way to evaluate A^J^^a2, contributing to 



(10) 



Aaa = - 



sin Qw 



[Ren'3T(s)-n''''^(0)] 
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which has been proposed long ago in [12]. The surprising fact is that the evaluation of 
a2 does not require to separate all individual flavor contributions to recombine them in 
the proper way. In fact, up to pcrturbative or very small contributions the hadronic shift 
of a2 is proportional to the self-energy correlation amplitude 11'^''' where 3 refers to the 
3rd component of the weak isospin current and 7 to the electromagnetic current. For 
the non-perturbative low energy range, it implies that the contribution corresponding 
to the u, d and s flavors actually requires no flavor separation in the SU(3) limit. This 
makes it possible to calculate Aa2 reliably, because the other heavier flavors may be 
safely separated by relying on pQCD weighting. The assumption is that the for Nf > 3 
the Nf — 1 lighter flavors above the Nf flavor threshold can be evaluated by pQCD. A 
detailed discussion of the approximations made is given in Appendix C of [12]. Given 
Hyj^ = n^J^^) + n^J) + n^^] for the continuum and n7^'^, ~ H^ + H" + H"^ + H"'/^ + H^ 
for the narrow resonances, we have the relations 

(11) n^;^,, c^lur, , + - ur, + - n/T; 

V / con 2 («rfs) ' g (c) ' 4 (b) 

for the background contribution and 

(12) n^i^in^ + ^n^ + ln^/^ + ^n^ 

\ / res 2 4 8 4 

for the resonance contributions. The p — uj mixing contribution usually included in the 
n'' taking into account the isospin / = component cj — > tttt in the 7 — > tttt — > 7 
channel is to be subtracted via the Gounaris-Sakurai parametrization (by setting to zero 
corresponding the mixing parameter). The coupling 02 can be "measured" in a charged 
current channel via Myy [g = 32): 

or via the neutral current channel sin Qf. In fact here running sin Qf{E) connects the 
LEP scale mixing parameter to the one of low energy v^s scattering 

(14) Sin^ Qe{Mz) = < . . , + A^^e,vertex+box + AKe,vortcx \ siu^ Qu^e ■ 

The first correction from the running coupling ratio is largely compensated by the 
z/p charge radius which dominates the second term. The ratio sin Q^^^/sin Og is 
close to 1.002, independent of top and Higgs mass. Note that errors in the ratio 
(1 — Aa2)/(1 — Aa) can be taken to be 100% correlated and thus largely cancel. 

Above result allow us to calculate non-perturbative hadronic correction in 77, jZ, 
ZZ and WW self energies. Gauge boson self-energies potentially are very sensitive to 
new physics (oblique corrections), which, however, may be obscured by uncertainties 
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Fig. 2. - sin^ Qw{Q) as a function of Q in the space-like region. Hadronic uncertainties are 
included but barely visible. Uncertainties from the input parameter sin^ 6w{0) ~ 0.23822(100) 
or sin^ 6w{Mz) = 0.23156(21) are not shown. Future ILC measurements at 1 TeV would be 
sensitive to Z' , H^^ etc. 



of the non-pcrturbative hadronic effects. For complete analytic expressions for clec- 
troweak parameter shifts at one-loop see [13, 14]. Another interesting version of running 
sin Qw{Q'^) one finds in polarized Moeller scattering asymmetries as advocated by Czar- 
necki & Marciano [15] (see also [16]). It includes specific bosonic contribution AKb{Q^) 
such that 



(15) 



f^is = -Q2) 



1 - Aa{s) 



AKb(g2) - AKb(O) 



where(*) , in our low energy scheme, we require k{Q^) = 1 at Q^ = 0. Explicitly [15], 



(16) AKbiQ^) 



2ns 



w 



in Cw H In t — 1 

12 18 V2 ^ 



(7 - 4z) cw 



+ -(l + 4z) 



z j r In^ + z (2 — z) In" ^ 



n7^ A ^n^ " f 42c^v + l 1 6cw + 7 



(*) Here Aa =dggvap(s,0.dO) and Aa2 =degvap(s,0.dO) are provided by functions from the 
package alphaQED. 
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with z = M'^/Q'^, r = a/1 + 4z, ^ = ^^, sw = sin^ Ovy and c^i/ = cos^ &w Re- 
sults obtained in [15] based on one- loop perturbation theory using light quark masses 
ruu = rud ~ nig ~ 100 McV arc compared with results obtained in our non-perturbative 
approach in Fig. 2. 

4. Adler function controlled split data vs pQCD 

A strategy to exploit the rather precise perturbative QCD predictions in a optimal 
well controlled way is to monitor QCD predictions via the Adler function D{Q^) in the 
Euclidean region by comparing theory and data there. 

(18) Di-s) = ^ s ^A«had(5) = - (12^2) s ^^ , 



OO 



(19) D{Q^) = Q2 f ^^ 



Ris) 



Low energies, resonances and thresholds prevent us from making reliable and precise 
predictions of R{s) in pQCD. Locally deviations between data and i?-prcdictions can be 
huge. In contrast, the smooth function D{Q^) is easy to compare and deviations show up 
at low energies only. A detailed inspection of the time-like approach shows that pQCD 
works well in "perturbative windows" like 3.00 GeV - 3.73 GeV, 5.00 GeV - 10.52 GeV 
and 11.50 GeV - oo. In the space-like approach pQCD works well for \/Q^ = —q'^ > 2.5 
GeV [17, 18]. Theory is based on results by Chetyrkin, Kiihn et al. [19, 20]. One thus 
requires data to calculate 

(20) Aahad(-so) = ^ r dQ''^^^ 

up to So = (2.5 GeV)^. Equivalently, Aahad(— so) can be directly calculated by (6) and 
used in (3). One obtains [17, 18] 

Aahad(-'So)'^''''' = 0.007337 ±0.000090 
Aahad(-^l) = 0.027460 ±0.000134 
A4ad(Af|) = 0.027498 ±0.000135 . 

The result includes a shift +0.000008 from the 5-loop contribution. The error ±0.000103 
in the perturbative part is added in quadrature. QCD parameters used are as{Mz) = 
0.1189(20), mdmc) = 1.286(13) [Mc = 1.666(17)] GeV, and mb{mc) ^ 4.164(25) [Mb = 
4.800(29)] GeV based on a complete 3-loop massive QCD analysis [21]. The latter 
results are in agreement with results from lattice QCD [22, 23, 24]. Results based on the 
Adler controlled split are A4^2iro„s(^l) = 0.027498 ± 0.000135 [0.027510 ± 0.000218] 
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Table I. - How much pQCD? AaJ^''jj(M|) x lO" pQCD part only. 



Method 


range [GeV] 


pQCD 






Standard approach: 
My choice 


5.2 - 9.5 
13.0 - oo 


33.50(0.02) 
115.69(0.04) 


— > 


149.19 (0.06) 


Standard approach: 
Davier et al. 


2.0 - 9.5 
11.5 - cxj 


72.09(0.07) 
123.24(0.05) 


— > 


195.33 (0.12) 


Adler function controlled: 


5.2 - 9.5 

13.0 - oo 

-oo - -2.5 

-Mz -^ Mz 


3.92(0.00) 

1.09(0.00) 

201.23(1.03) 

0.38(0.00) 


— >• 


206.62 (1.03) 



or a-\Ml) = 128.962 ± 0.018 [128.961 ± 0.030] in braces for comparison the results 
obtained by the standard approach. 

A comparison of error profiles between Aal^^^^M"^), Aal^^^{—SQ) and a^ may be found 
in [18]. Note that our approach, with a conservative cut of ^/sq = 2.5 GeV, does not 
rely substantially more on pQCD than standard analyses by Davier, Hocker et al. [9] and 
others (see Tab. I). Further progress is possible due to progress in methods to include 
the hadronic r-dccay data [25, 26]. 

5. The FORTRAN package alphaQED 

The FORTRAN package alphaQED . tar . gz [1] for calculating the SM effective cou- 
plings includes two versions: 



• 



alphaQEDreal [FUNCTION funalpqcd] providing the real part of the subtracted 
photon vacuum polarization including hadronic, leptonic and top quark contribu- 
tions as well as the weak part (relevant at ILC energies). Hadronic, leptonic. top 
and weak contributions are accessible separately via common blocks 

common /resu/dalept ,dahadr ,daltop,DalphaweaklMSb 
common /resg/dglept ,dghadr ,dgetop,Dalpha2weaklMSb 

alphaQEDcomplex [FUNCTION funalpqedc] provides in addition the correspond- 
ing imaginary parts. See Fig. 3. 

corresponding options alpha2SMreal and alpha2SMcomplex are available for the 
SU{2)l coupling a2 = .9^/47r. 

The functions are available for the space-like and the time-like region. The complex 
versions require to install the Rhad package of Harlander and Steinhauser [2] (FOR- 
TRAN package version rhad-1.01 (March 2009 issue)). The latter also provides the QCD 
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• 
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Qem(E) as a complex function 




INTERFACE 
rhad subdir V rhadHS VH^ 



Fig. 3. - Structure of alphaQEDcomplex. The corresponding diagram for alphaQEDreal is much 
simpler as it involves the upper part only. 



coupling q;3(s) = as{s)- The imaginary part given by the bare R^^'{s) is provided in 
parametrized form by Chebyshev polynomial fits. For sample plots I refer to the package 
description on my web page http://www-com.pliysik.hu-berlin.de/~fjeger/. The 
"organigram" of the program is shown in Fig. 3. 
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